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INTRODUCTION 


The  development  of  durable  damage  tolerant  composite  air¬ 
craft  structures  requires  a  sound  understanding  of  the  effects 
of  defects.  These  effects  are  qualitatively,  as  well  as  quanti¬ 
tatively,  different  from  the  fracture  mechanics  of  similar  metal¬ 
lic  structures.  In  recent  years  various  studies  have  been  con¬ 
ducted  for  assessing  the  effects  of  various  kinds  of  defects  in 
composite  components.  These  studies  have  reached  a  point  where 
an  effort  to  integrate  various  advances  is  now  desirable  in  order 
to  define  a  reliable  useful  methodology  for  defect  evaluations. 

This  capability  is  needed  for  fleet  life  management  and  repair/ 
replace  maintenance  decisions.  To  meet  this  objective  the  fol¬ 
lowing  major  tasks  were  defined  for  the  present  program. 

I.  Formulation  of  a  Master  Plan  -  Elements  of  the  methodology 
for  generalized  fracture  mechanics  and  NDE  of  defect  criticality 

in  practical  aircraft  structures  are  to  be  studied.  The  study 
should  address  not  only  such  technical  issues  as  type  of  defects 
and  method  of  evaluation,  but  also  such  management  issues  as  tech¬ 
nology  transfer  and  training.  Technology  gaps  in  various  areas 
should  be  identified  and  a  program  plan  to  fill  these  gaps  should 
be  drawn . 

II.  Perform  analytical  studies  and  experimental  investigation 
for  certain  types  of  defects  which  are  extremely  critical  like  de¬ 
laminations,  bonded  joints,  etc. 

This  report  summarizes  the  yorks  performed  to  date  in  various 
areas  during  the  first  year  of  the  present  contract. 


PRELIMINARY  MASTER  PLAN 


An  attempt  is  made  in  this  plan  to  identify  various  aspects  of 
the  operational  problem  of  nondestructive  evaluation  of  criticality 
of  defects  in  graphite  epoxy  composite  components.  Available 
nondestructive  test  methods  and  criticality  assessment  techniques 
are  evaluated  and  a  program  plan  describing  the  tasks,  which 
should  be  performed  for  effective  application  of  NDE  to  composite 
parts,  is  described.  The  plan  was  submitted  in  March,  1983,  and 
was  presented  in  the  symposium  on  NDE  of  Criticality  of  Defects 
in  Composite  Laminates  held  in  Valley  Forge,  PA,  on  May  23-24, 

1983,  organized  by  Materials  Sciences  Corporation  under  the  spon¬ 
sorship  of  Naval  Air  Development  Center  as  a  part  of  the  present 
program.  Based  on  the  discussions  and  other  presentations  at  the 
meeting  as  well  as  further  review,  the  plan  was  revised  and  re¬ 
submitted  (see  ref.  1) .  Work  on  an  expanded  plan  identifying  a 
detailed  program  to  fill  up  some  of  the  technology  gaps  is  being 
continued. 
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SYMPOSIUM  ON  NDE  OF  DEFECT  CRITICALITY 


As  described  in  the  last  section  a  limited  attendance  sym¬ 
posium  was  conducted  to  discuss  usefulness  and  shortcomings  of 
the  NDE  methodologies,  which  can  be  currently  used  in  practice, 
areas  where  further  research  and  developmental  work  is  needed 
as  well  as  the  present  state  of  the  art  of  inspection  and  repair 
technology  of  composite  components.  Vu-graphs  of  presentations 
and  transcribed  discussions  are  being  compiled  in  the  form  of  an 
informal  symposium  proceedings.  This  work  will  be  completed  by 
October,  1983.  A  discussion  for  some  of  the  presentations  is 
given  next. 

1.  "Use  of  Thermography  to  Detect  and  Quantity  Damage  in 
Fiber  Composite  Laminates",  by  P.  V.  McLaughlin,  Jr.  Monitoring 
the  surface  temperature  due  to  an  externally  applied  thermal 
field  by  the  use  of  an  infrared  camera  is  demonstrated  to  be  a 
useful  tool  in  locating  and  quantifying  defects.  From  the 
point  of  view  of  practical  application  there  are,  however,  some 
problem  areas  like  heat  generation  and  noise.  Increased  sen¬ 
sitivity  of  the  monitoring  instrument  may  be  helpful  to  sort 
out  some  of  the  difficulties. 

2.  "Ultrasonic  Correlation  to  Porosity  in  Graphite  Epoxy 
Laminates",  by  R.  Collins.  This  subject  is  of  great  importance, 
since  quantifying  porosity  by  NDE  is  a  difficult  task.  Although 
the  results  indicate  that  some  correlations  can  be  established 
between  porosity  as  determined  from  image  analysis  and  the  Db 
gain  at  which  an  indication  fails  to  trigger  the  recorder  am¬ 
plitude  gate,  further  work  appears  necessary  before  this  can  be 
used  in  routine  applications. 

3.  "Hole  Wear  in  Composites  due  to  Bearing  Loads",  by  R. 
Badaliance .  Fastener  hole  wear  damage  due  to  fatigue  loading  is 
found  to  be  a  local  one  caused  by  matrix  cracks  and  subsequent 
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pulverization  of  the  material  due  to  bearing  loads .  A  semi-em— 
pirical  correlation  between  strain  energy  density  stored  in  an 
infinitesimal  volume  element  at  the  location  of  maximum  bear¬ 
ing  pressure  and  constant  amplitude  hole  wear  lives  is  attempt¬ 
ed.  This  relationship  and  Miner's  rule  are  utilized  to  compare 
predicted  spectrum  lives  with  test  data  for  different  laminates 
under  various  environments .  The  approach  appears  to  be  useful  as 
a  design  tool.  However,  its  usefulness  as  a  NDE  technique  is 
limited,  since  no  correlation  is  attempted  between  extent  of 
damage  and  residual  life  time. 

4.  "On  the  Role  of  Damage  Tolerance  Analysis  in  Nondestruc¬ 
tive  Evaluation  of  Composite  Structures",  by  D.  J.  Wilkins.  Ap¬ 
plication  of  semi-empirical  growth  laws  and  3-D  finite  element 
analysis  for  calculation  of  energy  release  rates  in  predicting 
growth  of  a  delamination  under  fatigue  loading  is  illustrated. 
Advantages  of  fabric  and  3-D  reinforcements  for  increasing  damage 
tolerance  of  composite  structures  are  also  pointed  out. 

5.  "Ultrasonic  Waves  and  Composite  Material  Testing",  by 
W.  Sachse.  Use  of  ultrasonic  attenuation  and  velocity  measure¬ 
ments  in  detecting  material  inhomogeneities  and  flaws  is  examined. 
Further  research  is  needed  before  these  techniques  can  be  used  in 
practice. 

6.  "Computer  Aided  Ultrasonic  Measurement  of  Tangential 
Flaw  Boundary  Growth  in  Composite  Materials",  by  R.  A.  Blake,  Jr. 

The  proceaure  of  obtaining  digital  ultrasonic  C-Scans  on  a  color 
graphic  computer  is  described  and  its  usefulness  in  monitoring 
growth  of  delaminations  and  other  kinds  of  flaws  is  discussed. 

The  technique  of  computer  aided  image  enhancing  appears  attrac¬ 
tive,  but  more  work  is  needed  for  routine  application.  It  will 
be  of  interest  to  compare  such  measurements  with  those  from  other 
techniques  like  penetrant  enhanced  X-ray  radiography. 

7.  "Interlaminar  Stress  Field  as  Influenced  by  two  interacting 
flaws",  by  A.  S.  D.  Wang.  Finite  element  calculation  of  stress 
fields  near  interacting  intralaminar  cracks  and  interlaminar  de- 
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laminations  and  energy  release  rates  at  such  crack  and  delamination 
fronts  are  reported.  Correlation  of  predictions  obtained  from 
such  analyses  and  test  data  appears  necessary  for  demonstrating 
their  usefulness. 

8.  "Repair  Technology  for  Advanced  Composites",  by  E.  Rosen- 
zweig.  Currently  used  NDE  procedures  and  repair  techniques  are 
discussed.  It  is  reported  that  depot  level  inspections  are  being 
de-emphasized,  placing  more  importance  on  on-site  inspection. 

This  might  necessitate  a  re-examination  of  currently  used  NDE  tech¬ 
niques  . 

9.  "Presentation  of  a  Preliminary  Plan  to  fill  the  Technology 
Gaps",  by  S.  N.  Chatter jee.  The  preliminary  plan  discussed  in  the 
previous  section  is  presented. 

10.  "Potential  of  Accoustic  Emission  Technique  for  Monitoring 
Damage  During  Quasi-static  and  Fatigue  Loading"  by  J.  Awerbuch. 

Use  of  this  technique  as  a  research  tool  is  demonstrated.  Its  use¬ 
fulness  as  a  NDE  technique  for  evaluating  criticality  of  defects 
appears  limited. 

11.  “Defect  Enhanced  Fatigue  Failure  in  Graphite/Epoxy  Com¬ 
posites",  by  R.  J.  Richards-Frandsen,  et,  al.  Delamination  induced 
failure  is  studied  experimentally  and  a  semi-empirical  correlation 
is  attempted.  Use  of  more  refined  analytical  procedures  in  cal¬ 
culation  of  stress  intensity  factors  may  yield  better  correla¬ 
tions. 

12.  "Effects  of  Porosity,  Delamination  and  Low  Velocity 
Impact  Damage  on  the  Compressive  Behavior  of  Graphite/Epoxy  Lam¬ 
inates",  by  R.  L.  Ramkumar.  Test  data  and  some  analytical  re¬ 
sults  for  different  types  of  defects  are  surveyed  and  their  use¬ 
fulness  for  design  applications  are  pointed  out.  Severity  of 
different  kinds  of  defects  is  also  discussed. 

13.  "Compressive  Strength  of  Composite  Laminates  with  Inter¬ 
laminar  Defects"  by  J.  W.  Gillespie,  Jr.  et .  al.  Growth  of  de¬ 
laminations  in  presence  of  imperfections  of  wavy  shape  under  com¬ 
pressive  loads  is  studied. 
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14,  "Damage  Zone  Development  in  Notched  Graphite  Epoxy  Com¬ 
posites"  by  J.  Slepetz.  The  complicated  nature  of  damage  growth  in 
notched  composites  is  discussed. 

In  addition  to  the  presentations  discussed  above  use  of  ap¬ 
proximate  methods  of  stress  analysis  were  considered  by  D.  Y. 

Konishi  and  L.  W.  Rehfield.  Overviews  of  work  sponsored  by  Air 
Porce  were  presented  by  G,  Sendecky j ^  who  described  some  attempts 
of  using  CAT-Scans  for  detection  of  defects.  D.  Mulville  provid¬ 
ed  an  overview  of  Navy  sponsored  works  in  different  areas.  From 
the  presentations  and  discussions  it  appears  that  NDT  methods 
and  criticality  assessment  techniques  (in  many  cases  semi -empirical 
ones  based  on  test  data  on  similar  materials)  are  available  for  a 
variety  of  defects.  The  problem  of  choosing  the  fight  tool  or  tools 
for  effective  utilization  of  NDE,  however,  requires  further  research 
and  developmental  work  in  this  direction.' 


-6- 


ANALYTICAL  TASKS 


Analytical  tasks  to  be  performed  under  the  program  are  de-- 
scribed  below. 

1.  Stress  analysis  near  multiple  elliptic  disbonds  in  a 
laminated  plate  and  correlation  with  results  from  experimental 
tasks  on  growth  of  such  delamina j:ions  under  cyclic  loading. 

2.  Analysis  of  delaminatioii  near  ply  drops. 

3.  An  analysis  of  failure  of  bonded  joints  and  correlation 
with  test  data  from  experimental  program. 

4.  Comparison  of  results  from  two  analyses  of  the  delam¬ 
ination  problem. 

5.  Demonstration  of  usefulness  of  analysis  methods  with 
application  to  a  practical  situation. 

A  method  of  analysis  has  been  formulated  for  the  multiple 
delamination  problem.  A  cluster  of  delaminations  are  considered 
as  shown  in  figure  1.  Details  of  the  analysis  are  given  in  ap¬ 
pendix  A-1.  A  computer  program  is  being  developed.  Some  progress 
has  been  made  in  the  bonded  joint  analysis  considering  possible 
separation  near  the  ends  in  peel  mode  but  allowing  for  shear 
transfer  through  an  interdigitated  interface.  Details  will  be 
reported  in  the  next  reporting  period. 
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EXPERIMENTAL  PROGRAM 


The  experimental  study  has  been  separated  into  three  test 
programs  with  five  individually  identifiable  tasks,  depicted  in 
Table  1.  Test  series  involving  multiple  delaminations  in  Mode 
II  shear  are  similar  to  those  conducted  under  previous  programs, 
(ref.  2,  3)  while  the  investigation  of  biaxial  flaw  criticality 
and  the  analysis  of  bonded  composite  structures  with  adhesive 
delaminations  are  new  areas  of  interest. 


TASK  li  -  MULTIPLE  DELAMINATIONS  IN  MODE  II  SHEAR 

The  effects  of  multiple  delaminations  on  Mode  II  dominated 
flaw  growth  are  under  investigation  using  flexural  fatigue  speci¬ 
mens  in  three-point  bending.  Experimental  procedures  have  been 
developed  in  previous  test  programs  investigating  Mode  II  flaw 
growth  for  single  elliptical  delaminations  (ref.  2,  3) .  The  three- 
point  bend  test  fixture  employed  is  shown  in  Figure  2,  where  de¬ 
laminations  are  implanted  at  equal  distance  from  end  support  and 
center  load  on  both  the  left  and  right  halves  of  the  test  specimen. 

A  total  of  30  specimens  are  being  tested  in  three-point  bend¬ 
ing  to  examine  the  Mode  II  dominated  response  of  delaminated  plates 
Twenty-four  samples  are  under  test  in  flexural  fatigue  and  six  spec 
imens  were  failed  in  a  static  mode  to  provide  baseline  data.  Table 
2  describes  the  sample  designations  and  test  specimen  configuration 
for  Task  1. 

Two  elliptical  debond  geometries  are  investigated,  b/a  =  1.0 
and  1.5,  where  the  minor  axis  "a"  is  1  inch,  and  normal  to  the  long 
itudinal  direction  of  the  specimen.  The  circular  debonds  have  been 
chosen  for  investigation  for  direct  correlation  with  data  gener¬ 
ated  previously  for  the  single  circular  debond.  Delaminations  are 
implanted  at  the  midsurface  (Location  1)  and  at  quarter- thickness 
in  both  the  compressive  (Location  2)  and  tensile  regions  (Location 
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Table  1 


Program  Outline 

a)  3  Test  Programs 

Multiple  Delaminations  in  Mode  II  Shear 

Biaxial  Testing 

Bonded  Composite  Structures 

b)  5  Identifiable  Tasks 

1 )  Multiple  Delaminations 

2)  Off-axis  Shear  -  Biaxial 

3)  Honeycomb  Crossbeams  -  Biaxial 

4)  Bonded  Joints  with  Defects 

5)  Bonded  Joints  without  Defects 
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Table  2 


Task  1 

Laminate  Stacking  Sequences 


A:  (0g/+45^/-45^/05/+45^/ -454/02)^ 


B=  (0^/1452/4452/0^)2^ 


C:  (-45^/+454/02o/-454/+454/Og)5 

Sample  Designations 


Ex :  1 .0-A  1 


1.2-B  3 


Task  "1" 

Defect  Configuration  "0",  none 
Laminate  "A" 

Sample  "1" 

Task  "I” 

Defect  Configuration  "2" 
Laminate  "B” 

Sample  "3" 


Multiple  Delamination  Configurations 


• 

Location* 

Configuration 

Left 

Right 

1 

b/ a 

=  1 

.5 

Location 

1 

Location 

1 

Location 

2 

Location 

3 

2 

b/a 

=  1 

.0 

Location 

1 

Location 

1 

Location 

2 

Location 

3 

3 

b/a 

=  1 

.0 

Location 

1 

Location 

1 

b/a 

=  1 

.5 

Location 

2 

Location 

3 

4 

b/a 

=  1 

.5 

Location 

2 

Location 

2 

Location 

3 

Location 

3 

*Location 

1 : 

Debond 

at 

midsurface 

compression 

Location 

2: 

Debond 

at 

quarter-thickness  in 

Location 

3: 

Debond 

at 

quarter-thickness  in 

tension 
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3)  of  the  flexure  specimen.  Fatigue  tests  are  run  at  a  single 
S-level  of  0.5  and  three  replicates  are  being  tested  for  each 
geometry. 

Table  3  presents  static  test  results  for  the  no-defect  con¬ 
figuration.  These  data  were  used  to  set  the  S-level  for  all  sub¬ 
sequent  fatigue  tests. 

For  both  the  A  and  B  laminates,  the  failure  is  initiated  at 
the  compressive  surface.  However,  the  fracture  patterns  of  the 
two  stacking  sequences  are  different  (see  Figure  3) .  The  B  lam¬ 
inate,  with  +  45^  plies /  shows  only  a  compressive  failure  of  the 
outermost  0°  plies.  The  A  laminate,  however,  has  both  a  compres¬ 
sive  failure,  as  well  as,  considerable  interply  fracture  between 
the  +  45^  and  -  45^  lamina.  In  the  B  laminates  (+  sub¬ 

laminate  appears  to  have  arrested  this  failure  mode.  Fracto- 
graphic  analysis  of  these  failure  surfaces  remains  to  be  perform¬ 
ed  in  the  second  year.  It  should  be  pointed  out  that  the  failures 
are  in  flexural  mode  and  not  induced  by  flaw  growth  because  of 
small  delamination  sizes.  Flaw  growths  are  noticed  under  cyclic 
loading  as  discussed  in  the  follov?ing  paragraph.  Large  size  dis¬ 
bonds  can  initiate  failure  due  to  catastrophic  flaw  growth  (ref.  4). 

The  fatigue  program  for  Task  1  is  progressing.  All  samples 
in  both  A  and  B  geometries  have  been  fabricated  and  C-scanned,  and 
all  of  the  A  laminate  specimens  are  being  actively  tested.  Samples 
from  geometry  1.1— An  have  been  tested  to  failure,  while  sample  sets 
1.2,  1.3,  and  1.4  are  underway,  all  having  undergone  at  least 
100,000  cycles,  with  many  exceeding  150,000  to  date. 

Ultrasonic  C-scanning  is  being  used  to  nondestructively  mon¬ 
itor  the  growth  of  the  delaminations  during  load  cycling.  In 
order  to  individually  monitor  the  flaw  growth  of  the  disbonds  im¬ 
planted  at  the  quarter-planes  and  the  midsurface,  four  separate 
C-scans  are  necessary  for  1.1,  1.2,  and  1.3,  specimens.  These 
scans  are  designated  as  C,  T,  MC,  and  MT,  for  compressive  quarter- 
plane,  tensile  quarter-plane,  compressive-viewed  midplane,  and 
tensile-viewed  midplane,  respectively.  Figure  4  describes  the  gating 
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Table  3 


Task  1.0  Results 


Max  Load 
(kg) 

Thickness 

(inin) 

Width 

(nun) 

1.0-A  1 

2800 

8.9 

76.6 

1.0-A  2 

2650 

•  8.6 

76.2 

1 .0-A  3 

2900 

9.0 

76.0 

Average 

2783 

1 .0-B  1 

2950 

9.0 

76.3 

1 .0-B  2 

2400 

8.1 

76.3 

1 .0-B  3 

2700 

8.7 

76.7 

Average 

2683 
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Figure  3 


Task  1  Static  Failures 
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and  viewing  scheme  employed  to  monitor  the  individual  growth 
patterns  of  the  delaminations,  since  disbonds  are  implanted  only 
at  the  quarter-planes  for  samples  1.4,  only  c  and  T  scans  are 
needed.  A  series  of  c-scans  for  one  sample  from  1.2,  1.3,  and 
1.4  configurations  comprise  Appendix  A2. 


In  earlier  static  tests  (see  figure  28  in  reference  4)  it 
has  been  noted  that  interlaminar  shear  cracks  "stepped  up"  through 
he  0  plies  in  which  they  were  implanted  to  run  along  a  0/45 
interface.  This,  however,  appeared  to  have  occured  after  the  dis¬ 
bonds  have  grown  sufficiently  large  (elongated)  in  the  original 
plane.  m  the  current  work,  disbonds  are  implanted  at  the  quarter- 
plane,  which  is  a  -45/0  interface  for  the  A  laminate.  To  examine 
the  behavior  of  the  crack  front  originating  between  -45/0  plies,  a 
series  of  C-scans  were  done  on  specimen  1.4-A  2  after  110,000 
cycles  using  a  narrow  gate  focussed  above,  at,  and  below  the  im¬ 
plant.  The  sequence  of  five  c-scans  are  presented  as  Figure  5 
The  ultrasonic  images  show  that  the  delamination  crack  front  does 
change  laminate  planes,  to  "step  up",  in  this  case  to  the  +45/-45 
interface,  cutting  through  -45°  plies. 

Flaw  growth  will  be  measured  radially  from  the  center  of  the 
implanted  defects  by  either  manual  or  computer-aided  inspection 
techniques  following  completion  of  the  test  program.  The  radius, 
r,  and  angle,  e,  measurements  are  expected  to  produce  superior  flaw 
characterization  as  compared  to  point  measurements  made  along  the 
axis  of  the  beam  as  has  been  performed  in  previous  studies. 


TASK  II.  -  BIAXIAL  TESTING  THROUGH  OFF-AXIS  SHEAR 

A  three  point  bend  test  using  an  off-axis  laminate  is  being 
used  to  generate  a  biaxial  state  of  Interlaminar  shear  stress.  The 
relative  magnitude  of  these  stresses  is  limited  and  dependent  upon 
the  material  properties,  layup,  and  the  angle  of  orientation  be¬ 
tween  the  loading  direction  and  principal  material  direction.  The 
effects  of  bending-twisting  coupling,  however,  are  known  to  die  out 
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with  increasing  number  of  plies.  The  specimen  geometry  and  test 
methodology  are  indentical  to  those  employed  in  Task  1.  The 
initial  stacking  sequence  -  laminate  A  -  is  a  15°  rotated  ver¬ 
sion  of  laminate  B  from  Task  1,  ie .  [+15^/ (+60/-30) 2] •  De¬ 
pending  on  results  of  this  orientation,  a  second  orientation  will 
be  defined,  most  likely  either  a  +30°  or  +10°  series.  Table  4 
shows  the  sample  configurations. 

In  all,  eighteen  static  and'  twelve  fatigue  specimens  are  to 
be  tested  in  off-axis  shear  to  examine  biaxial  loading  response 
of  graphite-epoxy  laminates  containing  delaminations. 

To  date,  static  specimens  from  laminate  "A"  have  been  tested 
to  failure,  and  the  fatigue  samples  have  been  fabricated.  Static 
test  results  are  shown  in  Table  5. 

The  dominant  failure  mode  for  the  static  test  specimens  is 
intra-ply  splitting  of  the  surface  plies  on  both  the  tensile  and 
compressive  sides.  Figure  6  shows  the  common  failure  topography. 
The  large  magnitude  of  surface  cracking  provided  many  pathways  for 
the  water  used  as  the  couplant  during  ultrasonic  inspection  to 
enter  the  laminate  and  mask  the  interlaminar  features.  To  cir¬ 
cumvent  this  problem,  the  laminates  were  dipped  in  common  para¬ 
ffin  wax  to  seal  the  surface  cracks.  The  excess  wax  was  scraped 
from  the  plates  and  the  samples  were  subsequently  scanned  in  the 
normal  manner.  A  C-scan  of  a  wax-sealed  specimen  is  depicted  in 
Figure  7,  which  shows  no  growth  of  disbonds  has  occured.  It  is 
expected  that  they  will  grow  under  cyclic  loading. 

TASK  III.  -  BIAXIAL  TESTING  USING  CROSS-BEAM  HONEYCOMB  SANDVJICH 
SPECIMENS 


A  simple  method  to  develop  biaxial  in-plane  stress  is  to  use 
a  cross-beam  specimen  as  shown  in  Figure  8.  The  material  under 
test  is  bonded  to  one  surface  of  a  high  density  honeycomb  core 
material,  a  stronger  skin  is  then  bonded  to  the  opposing  face, 
creating  a  sandwich  structure  in  the  form  of  a  cross.  The  per- 
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Table  4 


Task  2  -  Off-axis  Shear  Specimen  Configuration 


Laminate  A: 

Laminate  B: 


Samples : 


(+154/(+60/-30)2)34 

(+304/(+75/-15)2)s4 

2.0-A  1,2,3 
2.1 -A  1,2,3 


Defect  Configurations; 

0  -  no  defect 

1  -  circular,  r  =  1.8” 

2  -  elliptical,  a  =  1.7",  b  =  2.0" 
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Table  5 


Static  Test  Results  -  Biaxial  Thick  Beams 


Max  Load 
(Xg) 


2.0-A  1  1790 

2.0-A  2  1800 

2.0-A  3  1840 

Average  1810 

2.1-A  1  1800 

2.1-A  2  1825 

2.1-A  3  1840 

Average  1822 


Thickness  Width 


(nun) 

(nun) 

8.9 

76.2 

8.7 

76.4 

8.9 

76.8 

9.0 

74.9 

9.1 

76.7 

8.9 

76.7 
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Composite  face 


Figure  8. 


Cross  sandwich  beam  test  specimen 


pendicular  beams  are  then  loaded  in  four  point  bending,  in  this 
case  producing  a  tensile  stress  in  one  direction  and  a  compressive 
component  in  the  other. 

A  test  program  employing  an  eight-ply  quasi-isotropic  skin 
has  been  developed  to  verify  the  viability  of  this  test  method  for 
study  of  delamination  growth  in  biaxially  loaded  panels.  A  total 
of  twelve  specimens  are  to  be  tested  as  described  in  Table  6. 

Teflon  disks  of  one  inch  diameter  will  be  implanted  for  the  six 
beams  with  delaminations . 

To  date,  the  test  fixture  has  been  designed  and  fabricated 
and  some  of  the  graphite— epoxy  skins  have  been  layed— up  and  cured. 

No  complete  specimens  have  been  finished  or  tested. 

TASK  IV.  -  BONDED  COMPOSITE  STRUCTURES 

The  strength  and  delamination  growth  in  bonded  composite 
structures  is  under  investigation  utilizing  the  single  lap  joint 
test  specimen  geometry.  The  lap  joint  consists  of  two  eight-ply 
quasi— isotropic  adherends  and  American  Cyanamid's  FM300M  adhesive 
used  commonly  on  several  aircraft  structures.  The  test  program  in¬ 
cludes  tension— tension  fatigue  and  static  tests  to  address  delam— 
ination  growth  and  strength  respectively. 

A  total  of  27  specimens  are  to  be  tested  in  the  bonded  joint 
program,  nine  in  static  mode  and  18  in  tensile  fatigue.  Two  teflon 
disbond  geometries  will  be  employed,  as  well  as,  a  no— defect  series. 
Figure  9  shows  the  delamination  geometries  to  be  tested. 

A  preliminary  study  was  undertaken  to  determine  the  pristine 
joint  strength  as  a  function  of  bond  length  and  to  generate  mechan¬ 
ical  property  data  for  the  laminate.  The  initial  set  of  specimens 
tested  showed  no  sensitivity  to  bond  length,  even  at  the  extremes 
of  10  mm  to  100  mm.  New  adhesive  was  acquired  to  retest  the  bond 
strength;  sample  fabrication  is  continuing,  including  dessication 

of  the  graphite-epoxy  laminates  prior  to  bonding. 

A  few  samples  were  bonded  with  Tedlar  film  implants  to  study 
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Table  6 


Honeycoir±>  Sandwich  Specimens 


Fatigue  Static 

(sfO.5) 

No  Disbond  3  3 

1"  Disbond  3  3 
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BONDED  JOINT  DISBONDS 


Figure  9  Bonded  Joint  Test  Specimens 


* 


the  ultrasonic  image  response  as  the  disbonds  were  moved  through 
the  thickness  of  the  adhesive  bond  line.  It  appears  that  Tedlar 
is  nearly  transparent  to  ultrasonic  NDE,  at  least  at  20  MHz. 
Following  this,  all  subsequent  implants  are  FEP  Teflon  film. 

Quasi-isotropic  mechanical  data  are  presented  in  Table  7. 

The  stress-strain  response  was  monitored  using  electrical  resis¬ 
tance  strain  gages  and  all  specimens  were  tested  in  a  screw-driv¬ 
en  test  machine.  Commercial  graphite-epoxy  Hercules  ASl/3501-6 
is  used  for  the  entire  test  program. 
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Table  7 


Quasi-isotropic  Mechanical  Properties 
AS1/350T^ 


Ell  =  49.84  GPa 
Vi2  =  0.32 

=  560  MPa 
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APPENDIX  -  A1 


3-D  ELASTICITY  SOLUTION  FOR  MULTIPLE  ELLIPTIC  DISBONDS  IN  A  LAM¬ 
INATED  PLATE. 


Statement  of  the  Problem 

A  cluster  of  elliptic  disbonds  in  a  laminated  plate  of  large 
size  subjected  to  either  equal  and  opposite  surface  tractions  on 
both  faces  of  each  disbond  or  self-equilibrating  tractions  on  top 
and  bottom  surfaces  of  the  laminate  over  a  finite  area  is  consid¬ 
ered.  It  is  assumed  that  the  edges  of  the  plate  are  far  away  from 
the  delaminations,  so  that  the  plate  may  be  considered  infinite  in 
X,  y  directions  (fig.  1) .  The  laminate  consists  of  a  finite  number 
of  layers  all  of  which  are  homogeneous  and  elastic.  A  layer  may  be 
a  unidirectionally  reinforced  fiber  composite  or  an  isotropic  mate¬ 
rial.  The  unidirectionally  reinforced  materials  are  assumed  to  be 
transversely  isotropic  in  planes  perpendicular  to  the  fiber  direc¬ 
tions  which  lie  in  planes  parallel  to  the  laminate  surfaces,  but  may 
be  orientated  at  an  arbitrary  angle. 


Solution  of  the  Equations  for  a  Layer 


The  equations  of  three  dimensional  elasticity  for  each  layer  of 
a  given  laminate  are  given  by 

‘^jk,k  ^  °  ^jk  ^  ^jk£m^£,m  (2) 

Using  Latin  suffixes  j ,  which  take  values  1,  2,  3  and  Greek  suffixes 
a, 3  which  take  values  1,  2  substitution  of  (1)  in  (2)  yields. 


^^jaJl3  ‘^j3£a^  ,a3  ^j3£3^£,33  ° 


(3) 


Note  that  C  .  =  C 

gkjlm  np 
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where  (n=j  if  j=k)  or  (n=9-j-k  if  jj^k) 

{p=£  if  £=m)  or  (p=9-il-m  if  £7^m) 

Taking  Fourier  transforms  in  directions 

'V  00  00  . 


—  00  — oo 


where  x  =  (x^,X2)  and  ^  vectors 


(4) 


(5) 


Inverse  of  (5)  is  given  by 

(xr^Cg)  =  J"  J"  (5^X3)  e  ^-•.^d^3^dC2 

—  00  —00 

Substitution  of  (6)  in  (3)  yields 

^  *\j 


General  solution  of  (7)  is 

«£(i,X3)  =  E  £^^q^^^^3}Bq(i) 

q=i 


(6) 


(7) 


(8) 


where 


“  ^^ja£3  ^j3£a^  ja£3  ^  ° 


(9) 


q  =  1,  6 

X  are  therefore  the  roots  of 

q 

Det.  [Cj3j^3Aq  -  (Cj^^3  +  ^ j3  ]  =  0  (10) 

A  's  are  either  purely  imaginary  or  complex  and  occur  in  com- 

q 

plex  conjugate  pairs.  Therefore,  let 
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Vs  "  '■  g  -  2,  3  (1. 

For  layers  of  unidirectionally  reinforced  composites  which 
are  transversely  isotropic  in  the  plane  perpendicular  to  the 
fibers  the  roots  can  be  expressed  in  simpler  form  by  utilizing 
the  following  transformations  (fig.  1). 


cos  9  -'sin6 


sin9  COS0 


cos9  -sin6 


sine  cosf 


Where  6  is  the  angle  between  the  fiber  direction  and  x^  axis. 
This  yields 

Introducing  displacement  variables  in  x-,  ^  x^  direction. 

Where  x^^  coincides  with  the  direction  of  fibers  so  that 


cos9  -sin0  0  u- 


sine  COS0  0  u. 


0  1  u. 


one  obtains  for  transform  variables  u. 


^  ixf"'' 


Where , 
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(16) 


q=l 


The  roots  are  given  by 

o  •'^2  ''2 

A  ‘  *  ^2  1 


(17) 


Where 


2 

K  = 

q 


1 

2 


-y^ 


4C32/C11 


f  q  =  1/  2 


K  -  (03^3^022  -  ^12  ~  2Ci2^55^/^11^55 


(18) 


1^3  ^44/^55 


C.  .  being  the  engineering  stiffnesses  in  contracted  notations, 
direction  1  conciding  with  the  fiber  direction  and  and 

denoting  the  axial  and  transverse  sheer  moduli  respectively. 
The  six  roots  X  are  therefore  given  by 


X  =  -  iXq 

q 


,q  =  1,  2,  3 


and 
where ,  X 


(19) 


For  commonly  used  fiber  composites  the  roots  are  all  distinct 
■  > 

and  A  in  (16)  can  be  expressed  as 
iiq 


Iq 


=  ’  ^2q  =  V2  '  %  =  -"'^q  '  q  '  1'  2 


where,  “  *^22'^^q^ 


(20) 


^13  '^23  ”^3'  ^33  ^^2 
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and 


(21) 


^£(q+3) 


Hence,  it  follows  that  equation  (16)  can  be  written  as 


^."(C  ,  x^)  =  -y  2 

^  ^  ^  q=l 


%  +  %  ®q+3 


(22) 


For  a  layer  deforming  symmetrically  (c=l)  and  antisymmetrically 
(c=2)  about  the  midplane  X-=0  the  following  solutions  are  ob- 
tained  by  taking  for  c-1  and  =  -B^  for  c=2 ; 


A  ^  e  c 

u,  =  I  B  A,  (x  )  ;  £  =  1,  2 

q  iiq  iq  3 


(23) 


^  .A  3  c>  11^ 

1^3  ^3q  ^2q  ^^3^ 

Where  A_,  =  iA-,  ,  A^  being  given  by  (20)  and  the  functions 

3q  3q  dq 

c 

F  ,  c  =  if  2  are  as  follows 
mq 


Iq 


Case  1  (C=l) 


cosh  igX^ 


Case  2  (C=2) 
sinh  igX^ 


(24) 


2q 


sinh  A  x_ 


cosh  X  x_ 

q  3 


Relationship  of  Boundary  Tractions  and  Displacements 


The  transformed  tractions  at  the  boundary  z  -  +  -j-  of  the 
layer  are  related  to  the  transformed  displacements  at  z  =  +  y  by 

the  following  relations. 
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T  (+4-)H 


a3i(+  2  ) 


h 


^32^^  2  ^ 


c 

h 

lOoo  (+  -5- 

L 


- 

1 

- 

c 

c 

c 

c 

’"11 

A 

^12 

,A 

^13 

O/ 

u£  (+  - 

^) 

2 

c 

sym. 

^22 

^23 

U2(+  ■ 

2  ^ 

c  > 

c 

1 

^^33 

.  i"3<+ 

^) 

2 

(25) 


c 


=  [K"]  U  \+  -^) 


where , 
c  c 


A  x;  x'  (A*  -  A*)  +  52  (A’Xjfj  - 


c 

A  k 


^  A  A 


12  -Ss  «2  <^2*2  -  *: 


C  C,  '  r  '  '  -k  '  '  *.  J.  c^2,  *  *.  £ 

^  *^13  =  Ss  ["  ^1  ^3^^  ■  ^2  ^3^^  "  ^2^  ^2  ^2  ^2^ 


^  ^"22  =  ^=^55  ^3 


(26) 


c  c 
A 


^^23  =  ^44  ^2  ■  ^2^2^2^  '  ^^3^  ^  "  '^2^^3] 


C  C  A  ^  ^  ^  * 

^  >'33  =  <=55  ^3  <^2  -  ''1' 

c 


and  A  =  X3  (^2^2^2  "  ~  ^2^  ^3 


and  fj,  j  =  1,  2,  3  are  given  by 

Case  1  (c=l) 

f .  =  tanh  A 
3  q  2 


Case  2  (c=2) 
coth  X 

q  2 
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The  quantities  ^2  etc.  in  equation  (27)  can  be  ex¬ 


pressed  as 


Cos6  +  Kz  sin9 ,  ^2  "  oj.i*u  .  ^2 


sine  +  U  cose 


cos^e  +  ^2  sin^e  +  cosO", 

^2^  =  ^2  sin^e  +  ^2  ~  2  ^1^2 


(27) 


A  = 


cos^e  +  sin^e  1+  j:2 
2“  /  ^2 


2  2  ^ 
sin  e  +  cos  e 


+  2  q?2  - 


-  l\sine  cose 


Transformed  tractions  in  ^2,  coordinate  systems  are  re¬ 

lated  to  the  corresponding  displacements  by 


[T  (+  -|*^  ]  =  [K]  [U  (+  \)  ] 


(28) 


where 


+  Ic 

11  22 


11 


-  + 


11  “  ^22 


cos2e  - 


k^2  sin2e 


c,  c^ 

^11  "  ^22 


12 


sin26  +  k^2  cos2e 


13  ”  ^13  "  ^23 


(29) 


c  >•  c 

^11  ^22 


c,  c^ 
^11  "  ^22 


22 


cos2e  +  k^2  ®iJ^2e 


23 


^2  sine  +  k22  <^26 


c> 

33  ""  ^33 
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Now  a  general  case  of  loading  on  the  layer  can  be  written  as  a 
sum  of  symmetric  and  antisymmetric  loadings.  Let  the  trans-^ 
formed  tractions  '^2'  at  z  =  +  ^  be  denoted  as  , 

where  t,  ,  t„ ,  and  t^  are  in  +ve  x,  ,  x_,  x-,  directions  and 
[u,  ,  u„,  iu-]  at  z  =  +  -7p  be  expressed  as  U,  .  Let  the  cor- 
responding  quantities  at  z  =  -  -j-  be  written  as  T2  and  £2 
respectively.  Then  one  can  write 


V  11  12  'V 


^1 

= 

[K] 

^1 

+  [K]  £2 

v 

21 

a. 

22  'v. 

T 

2 

= 

[K] 

^1 

+  [K]  £2 

where 

"  1 

2.  1 

2  1 

2 

(kii 

+  (k^2 

+  k^2) (^13 

^13^ 

11 

1 

2  1 

2 

[K] 

^^22 

+  ^22^  (^23 

+  k23) 

12 

[K] 


21 

[K] 


sym. 


<^11  ■ 

^11^  ^^12 

1 

2  1. 

^^12 

^12^  ^^22 

1 

2-  1 

J^13 

^13^  ^^23 

12  ^ 

(k33  +  k 

33^ 

2  1 

2 

^12^  ^"^13 

^13 

2  1 

2 

^22^  ^“^23 

^^23 

2  1 

2 

^23^  (-k33  + 

^33 

(31) 


22  11 
and  elements  of  [K]  are  obtained  from  these  of  [K]  by  changing 

the  signs  of  ^^23'^^31'  Expressions  for  and 

are  given  by  (29  ),  (26)  and  (27). 
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Satisfaction  of  Interface  Conditions 


th  t  h 

At  an  interface  I  between  (il-1)  and  i  layer  {l<£.<n) 
without  any  disbonds  the  displacement  vectors  of  two  adjacent 
layers  above  and  below  the  interface  must  be  identical  and  the 
corresponding  traction  vectors  must  add  up  to  zero.  Note  that 
1  and  n  indicate  the  top  and  bottom  surfaces  of  the  laminate 
where  tractions  are  either  zero  or  have  a  prescribed  value.  Let 


—Z  denote  the  displacement  vector  at  all  such  interfaces. 
At  interfaces  (j  =  1,  2..m) ,  containing  disbonds 


let 

'V£-rl 

U  ^  ^ 

-2  -Z . 

11* 

O;  * 

V. 

—3 

J 

J 

3 

and 

^Z  . 

=  y-n . 

- 

22* 

Oj* 

V. 

-3 

D 

3 

: 

(32) 


The  superscript  on  U  in  the  left  hand  sides  of  (32)  indicates 
the  layer  and  the  subscripts  1  and  2  stand  for  top  and  bottom  sur¬ 
faces  of  that  layer.  Vj  in  the  right  hand  side  of  the  equation  is 

the  transform  of  the  displacement  discontinuity  vector  at  inter- 

ag* 

K, 


'y>Z  .-1  '^Z  . 
face  Z .  (=  U  ^  -  U,^) . 

^  aB 


X  is  the  asymptotic  form  of  the  stiff- 

X'  j 

ness  matrix  [K„  ]  of  layer  £ . ,  i.e.  the  limiting  form  as  the  layer 
j  ^ 

thickness  h  approaches  infinity.  Note  that  if  one  writes  ^^2.  ,^2^ 
[^cos(})  ^sin^)] 


then, 


aB*  otB 

[K  ]  =  ?[F(^)] 


(33) 


where  the  elements  ^22'  ^12  ^22,  ^33  behave  in  the 

following  manner 

f^^(TT  +  <|>)  =  ^11 

and  the  elements  f^^  ^23  the  relationship. 
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(34b) 


* 

The  matrix  [K^  ]  in  (32)  is  given  by 


0  -1  11* 

3  3 


22* 


(35) 


Since  the  tractions  on  the  top  and  bottom  disbond  sur¬ 
faces  are  either  zero  or  equal  and  opposite  the  traction  vectors 
at  such  an  interface  5,^  acting  on  the  adjacent  layers  add  up  to 
zero  -as  in  the  case  of  interfaces  without  any  disbonds.  In  view 
of  (32)  these  conditions  can  therefore  by  written  as 


V, 


[K^] 
3nx3  n 


3 


Xn 


=  [B] 
3nx3m 


^1 

• 

0 

• 

• 

0 

l\,* 

V. 

-3 

m 

+ 

•  1 

• 

• 

• 

• 

9 

_ 1 

(36) 


The  matrix  is  the  global  stiffness  matrix  in  the  trans¬ 
formed  domain  and  is  of  the  following  form 


rrrG, 

LjC'-  ] 


11 

K 


12 

K, 


1  “2 
22  11 
(K^  +  K2) 


sym. 


O 

12 

K. 


0 


0 

0 


22  11  12 
(K2  +  K3)  K3 


22 

K 

n 


(37) 


The  matrix  B  is  of  size  (3n  x  3m)  and  the  elements  of  columns 

3^..-2,  3£  .-1,  and  3£  .  are  given  by, 

3  3  3 
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0 


Rows 

31  .-3^31. -1 

1  3 

- 

12 

0  11 
l[K^l  IK^l 

(33) 

1 

Rows 

3il.-2^3)lj 

11 

0 

22*  22 

0  11*  1 

[K  ]  [K„  ]  ^ 

)6j  y.. 

Rows 

3^1  .+1^31  .  +  3 

3  3 

3 

21 

0  22* 

0 

Of 

T  and  T  are  the  transforms  of  applied  traction  vectors  on  top 

and  bottom  surfaces  of  the  laminate.  It  should  be  noted  that  all 

elements  of  B  in  (38)  approach  zero  expotentially  as  /^2  ^  ^2 

1  2 

approach  infinity.  ^ 

Equation  (36)  can  be  formally  solved  to  express  as 


m 


=  I 
p=l 


fD.  ]  V  + 
ilp  -p 


[E 


'h 


T 

— n 


(39) 


The  transform  of  the  traction  vector  on  the  top  face  of  layer 
(j=l^  2.  .  ra)  (corresponding  to  bottom  surface  of  disbonded  inter¬ 
face  i.)  can,  therefore,  be  written  (with  the  help  of  (30),  (32), 

and  (39))  as 


M  . 

11 

^Si  . 

1,^  - 

T.  = 
-3 

1 21^ 

r\j-k 

m 

[C^P] 

Of  * 

T  .  = 
-3 

*  ■ 

=  [C 

+  E 
p=l 

V  + 

— ^T', 

ir 

12  (40) 

+  (Kj  1  U  3 

3 

a<0 

T.  (41) 

-3 
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where 


(42) 


*  11*  0  22* 

[C.]  =-[K  ]  [K  ]  [K  ] 

3  x,j  x,j  ^3-1 


11 


12 


(43) 


11* 

11 

22* 

^D 

•  ^D 

12 

^j+1 

11* 

3 

^  ^^'^j+1^  *^(D+1)P 

6 .  being  the  Kronecker  delta.  The  term  T?  given  by 
IP  ^  3 


DP 


'X/O 

T  .  = 
-J 

[F-iJ 

-1 

[Fjjl 

% 

T 

-2 

11 

i2 

J 

.P-'  * 

]  [E 

j  (^j+l)P 


,];  P'  =  1,2 


(44) 


The  elements  of  [C^]  behave  like  (33)  whereas  these  of  [C^^] 
approach  zero  expotentially  as  . 


V.,  the  transform  of  the  displacement  discontinuity  vector 
“1 

a  interface  containing  nuriber  of  elliptic  disbonds  can  be 
expressed  as 


1  ^3*  ix  r 

^3  =  2^  ^  ^A.  ''j  'S)  e 

q=i  Dq 


dA 


(45) 


The  inverse  transform  of  equation  (41)  is  therefore, 

*  I  CO  ,oo  — ix* 

Tj(x')  =  L  Tj  e  -  -  d?^d?2 


Qj 


1  ^-^OOOO  ^“1 

=  — 2  s  /  /  [C  (x)  e 

4tt  q=l 


q' 

7  . 
-D 


i (x  -  x' ) . ? 
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q’ 


1  ni  Qp  ^  ^  . 

+  E  I  f  f  [C^P]  V  e 

4/  p=i  q=i  ^ 


i (x  -  x' ) .? 


d?2  dA 


-oo  -oo  ^  “IX  •  E 

T.  e  —  — 

2tT  “"00  ""OO  ^ 


+  —  /  / 


dSjdCj 


j  =  1,  2  . . .m. 


Since,  the  tractions  on  disbonds  at  interface  are  either 

zero  or  prescribed,  equation  (46)  is  a  set  of  integral  equations 

CT  ^ 

for  evaluation  of  the  unknown  functions  V?  ;  q  =  1,  2,  .  .  .  , 

j  =  1,  2,  ...m.  For  obtaining  numerical  ..solution  of  (46)  the 
elliptic  disbonds  are  mapped  on  to  circular  regions  and  equation 
(46)  is  discretized. 
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Appendix-A2 

Representative  C-scans  of  Multiple  Delamination  Specimens 

Gated  at  C,  MC,  MT,  T  (see  fig.  4) 
for  Different  Load  Cycles 
See  Figure  4  and  Table  2 
for  Gating  Procedure  and  Specimen 
Details 
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Specimen  1.2 -A3,  C,  MC,  MT,  T 
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Specimen  1.3-Al  C,  MC,  MT,  T 


TV- 


650  0 


55000  65000  85000  110,000 


Specimen  1.4-A3,  T, 


